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We characterise the spontaneous emission time and direction from small numbers of dipole-coupled 
two-level atoms (2LAs) in the presence of incident fields. We show how to use adiabatic passage 
to admit population transfer between states in the one-quantum subspace for two and three 2LAs. 
Our method is a multi-atom generalisation of stimulated-Raman-adiabatic-passage (STIRAP) for 
a single multi-level atom. We use numerical results to justify an ansatz that enables us to give 
' r S-T . analytical expressions for the directional emission which depends on the incident fields. Our results 

j£z£ • admit a characterisation of the efficacy of population transfer in small numbers of dipole-coupled 

' 2LAs, and are applicable to proof-of-principle experiments involving dipole-coupled 2LAs. 
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C/5 ! I. INTRODUCTION 

00 . 

Atomic systems have proven to be valuable for experimental demonstrations of entanglement and violations of Bell 
inequalities [l| , quantum information protocols @, H, 0| , cooperative effects such as superradiance S, and quantum 
beats 0. A pair of electric dipole-coupled two-level atoms (2LAs) has been especially well studied [3] , as this is the 
fS ,■ simplest atomic system for demonstrating these effects, and three 2LAs is the smallest number of atoms that encode 
a decoherence-free qubit (two-level system) Common to all demonstrations of quantum beats, cooperative effects, 
entanglement, and quantum processing tasks is the need for external control. In most physical systems, external 
control requires incident fields, which can alter the characteristics of the system. Quantifying the effect of such fields 
on the emission properties of small collections of 2LAs is one purpose of this paper. 

We focus on systems in which the atomic coupling is described by a dipole-dipole interaction. We assume the 
2LAs are closely separated, and so permit only global control. Previous work has quantified global preparation, 
manipulation, and readout of collective states [9|, [l0(. Here, we show how to transfer population between entangled 
states in both two and three 2LAs using resonant counter-intuitive pulses similar to stimulated Raman adiabatic 
passage (STIRAP) [ll[. There are two reasons why it is surprising that this technique is applicable. First, there 
is no null eigenstate in the field Hamiltonian — normally required for STIRAP — and second, the 2LAs are spatially 
separated. Our results exploit the energy-level splitting due to the dipole-dipole interaction, allowing certain (fast- 
oscillating) terms to be neglected. This enables an eigenstate of the field Hamiltonian to be approximated by a 
dark-state. The approximate dark state is between two maximally- entangled states, so we call our method entanglcd- 
STIRAP. For two 2LAs, entangled-STIRAP is possible at large separations, and so we hope is amenable with present 
\ technology. 

We then examine the spontaneous emission characteristics in the presence of global control-fields incident on two 
and three 2LAs. Quantifying the emission characteristics serves two purposes. First, the decay timescale allows for the 
detrimental effect of the control-fields on the emission timescale to be characterised. Second, the emission direction 
provides a means with which to quantify the efficacy of the transfer. This gives a practical way to measure the fidelity 
of the population transfer. 

The paper is structured as follows. In Sec. UH we describe the physical system, and a particular unravelling of 
the Lindblad master equation. Then, in Sec. IIII1 we focus on two 2LAs. In Sec. IIII A[ we propose two methods for 
coherent rotation in the one-excitation subspace: collective two-photon Raman transitions and entangled-STIRAP. 
In Sec. IIII Bt we use numerical results to justify an ansatz that enables us to find an analytical expression for the 
directional emission matrix element in the presence of an incident field. We then turn our attention to three 2LAs. In 
Sec. lIV Al we show how to use entangled-STIRAP to rotate information encoded in the slowest decaying excited-states 
in three 2LAs. Then, in Sec. IIV B[ we give an indication of the effect of the incident fields on the emission pattern 
and time, and we show numerically how the emission pattern is affected by the incident fields. 
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II. PHYSICAL SYSTEM 



Our analysis is in accordance with standard quantum optical methods for electric dipole-coupled two-level atoms 
(2LAs): a quantum master equation employing the rotating- wave and Born-Markov approximations. We assume that 
the spatial extent of the 2LA is much less than the resonant wavelength A = 27r/fc , and so the 2LAs are effectively 
point dipoles with resonant frequency uj . We describe the driving field as a classical bichromatic field that drives all 
2LAs simultaneously, but due to the small interatomic separations cannot drive each 2LA individually. The raising 
operator for atom i is <7j_|_ = |l)i(0| = crt_ and &iz = [o-j+,<7j_]. The transition-matrix element of the atom is given 
by d = j(0| di with the dipole operator di for atom i. The dipoles of the atoms are identically oriented, with 

d ■ fij = cos a for d the unit vector in direction d, and fij = rtj jr^ the unit vector in the direction of the separation 
vector rij between atoms i and j. Since all the atoms are identical, the matrix elements of their dipole operators are 
equal: d = djVi. 

The free evolution of N dipole-coupled atoms, including all non-ncarcst-ncighbour interactions, without an incident 
field is [h = 1) 

N N N 

i— 1 i—1 



with 

3 7 e if 



[4 sin 2 a - (1 - (1-3 cos 2 a)] , (2) 



'j 



for i,j = 1, . . . , N, £ij = k$rij, and 7 the single-atom decay rate. By symmetry, S»j = 'Bji, and we define = 
Re{Sy } and 7^ = — 2 Im{S,j}. The coefficients Ay and 7y correspond to the dipole-dipole interaction and the actual 
process of photon emission respectively. The i th atom is located at ri = (i — 1)sXqZ, for z the unit vector along the 
interatomic axis, and s the interatomic separation in units of emission wavelength. 

The laser has a bichromatic electric field E(r) = E^ir^+Evir), with E^r) and E v (r) the electric field amplitudes. 
E(r) interacts with the atoms via a dipole coupling and so we introduce the Rabi frequencies — d- E^c~ lkti ' ri and 
£„ i = d- E y c~ lk " ri for wavevectors and k u . Within the rotating- wave approximation, the interaction Hamiltonian 
is 

N 

H I =J2£ i a i -+R.c, (3) 

i=l 

for which Si are described by a time-dependent bichromatic external field: 

Si S^^iG ^ ~\~ Su^iQ . (4) 

Due to the different positions of the atoms, the field Si differs for distinct atoms. The total effective Hamiltonian for 
the no-jump evolution is H c s = Hg + Hj. 
The jump operators are written [l3| 

N 

S{6,4>) = y/-yD{6, <f>)dfL ^ e -ik R(e,<f,yv j& ._^ ^ 
i=i 

which apply when a photon is detected in the far field within the element of solid angle dil in direction R(0, 4>). The 
dipole radiation pattern is 

D(M) = J-{l-[d-£(M)] 2 }, (6) 
which is the directional emission from an isolated atom. So, the master equation is written 

p = -i(H ctfP - P Hl s ) + J S(6, <f>)p&{8, 0), (7) 

which corresponds with the standard Lehmberg master equation for N atoms [HI, EH EH . Writing the master equation 
in this way is useful for analysing spontaneous emission in the presence of incident fields. 
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FIG. 1: Populations of the levels a,b,c,d when H a s from Eq. (|13p is applied to two atoms initially in |b) with (a) £ M = £ v = 37, 
= 30 7 , a = 0, and £12 = 1/5, and (b) £ M (t) = (3/4) 7 exp[-(l/5)(t - i M ) 2 /T 2 ], £„(t) = (3/4) 7 exp[-(l/5)(f - tuf /T 3 }, for 
T = 5.57- 1 , t M = 4.5T, i„ = t M + 2.75T, a = 0, w s = 0, and £12 = 1.25. 



When studying directional emission properties we only consider a = n/2. So, the jump operators can be written 



N 



S(9) = y/~/D(0) sin 8d8 e-' l2% ( j -V s cos °& j 
j'=i 



(8) 



where we have exploited the symmetry of the atomic system and integrated over azimuthal angle. In this instance, 
the dipole radiation pattern is written 



ore) = I - 1 sin 2 e . 

V ' 4 8 



(9) 



III. TWO ATOMS 

We define the basis of two atoms as |a) = 1 00) , corresponding with both 2LAs in the ground state, and |d) = |11) plus 
the symmetric state |c) = (l/\/2)(|10) + 1 01} ) and antisymmetric state |b) = (l/v / 2)(| 10) — 1 1) ) . The corresponding 
energies are 0,2wo,^o + A and ujq — A, respectively. The frequency difference between levels |a) and |b) is given by 
w a b, and similar notation is applied to other transitions. For two atoms, A = Rc{Si2}, T = — 2Im{Si2}, and the 
decay rates of the states |b) and |c) are 7 — T and 7 + T respectively. Note that the decay rate for the b-a transition 
approaches zero in the limit £ — ■> 0. The atoms are arranged according to r\ — and r 2 = r. 



A. Manipulation 



Bell state preparation of two dipole-coupled 2LAs has been studied both experimentally [17|,[T8| and theoretically [10|, 
IT9L [20l | . In particular, in Ref. [l(J | preparation of Bell states is performed by coherent manipulation via an incident 
laser field. We use the same preparation method here, and we ignore decay for the rest of Sec. IIII Al For preparation, 
only a single laser field is required. So, in the interaction picture with respect to the Hermitian part of Hs, the 
Hamiltonian H e s becomes 



H c s = 



1 

2Vf 



(^,1 + ^,2) (e- iAt |a> <c| + e iAt |c) (d| ) - (£ M - f„, a ) (e iAt |a) (b| - e - iAt |b) (d| )] ^~"o)t + H . c . 



(10) 



4 




FIG. 2: (a) Rabi frequencies of £ M = J2„ e s £ > where £ = 6O7 exp[-3.27(t - n) 2 ] and S = {1.23,3.31,4.53,6.62,...}, and 
= X^ngs' ^ f° r "S 1 ' = {1-72, 2.94, 5.02, 6.25, . . .} and (b) populations of the levels a,b,c,d when Rabi frequencies as depicted in 
(a) are applied to two atoms, initially in state |b) for a = 0, £ = 1/5, and ujs = 0. 



where £^1 and £^2 are the Rabi frequencies at atoms one and two respectively. To prepare |c), the incident fields 
are set to satisfy £^\ = £^2 = £p and the laser frequency tuned to lu^ = cj a b = loq + A, which gives 

i/ cff =^|(|a)<c|+e 2iA *|c)(d|)+H.c., (11) 

for which the coupling |c) (d| is rapidly oscillating for large A. This enables the Hamiltonian to be approximated by 

iTeff = ^= |a) (c| + H.c, (12) 

which allows high-fidelity oscillations between states |a) and |c). The antisymmetric state can be prepared using 
£^1 = — = £^. Note that this condition cannot be satisfied using a single laser field, because the distance 
between the atoms is much smaller than the wavelength Ao- One possibility is to use a standing- wave configuration, 
which is described in detail in Ref [h| . 

Once a single-excitation Bell-state is prepared, coherent Rabi oscillations within the single-excitation subspace can 
be performed using bichromatic fields. Here, we present two methods that admit this transformation, both of which 
rely on neglecting rapidly oscillating terms. We hope that these proposals will be of interest for proof-of-principlc 
experiments regarding coherent manipulation of atomic systems, particularly in the solid state. In particular, although 
the 2LAs are closely spaced, the proposal here is comfortably within the capabilities of experiments performed on 
nitrogen- vacancy centres in diamond [2l| . 

We use two fields that satisfy = £ Mj 2 = £^ and £ v .\ = —£ u ,2 = £v respectively. In Ref. [lCj| . an example of an 
experimental setup required to achieve these conditions is given. The field Hamiltonian in the collective basis can be 
written 

iTeff =-^(£ M e" itW5 |c) (dl+^e"^ |b) (d\+£^-^ 2A+ ^ |a) (c| - E^ 2 ^^ |a) (b| ) +H.c, (13) 

where uj^ = loq — A — ujg and o>„ = luq + A — tug for tog the detuning. This Hamiltonian is the same as that stated in 
Ref. [HI. 

Our first method applies the idea of single-atom two-photon Raman transitions to the collective system of two 
atoms. We show that coherent transfer between Bell-states is made possible by judicious choice of detunings and Rabi 



frequencies. Figure 1(a) shows the effect of the bichromatic fields taking into account all four levels. We adiabatically 



eliminate the off-resonant levels to give 

iJeff = £eff|c)(b|+H.C., (14) 
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FIG. 3: Average angular distribution of emitted photon with and without an incident field for a = 7r/2, £ = 1/5, and 
< t < 20y -1 , for emissions every jSt = 0.001, and for |^(0)) = |c). The distributions $b(0) and $<:(#) are stated in Eqs. pO)) 
and Q2ip respectively. The angular distribution including the incident fields <&n(8) is shown in Figure(a), for an exact numerical 
solution including all four levels (solid line) and for the ansatz proposed in Eq. (|22|l (dashed line) where Q, = tanhi. 



for the effective Rabi frequency 



£*£ 

£ " = it- (15) 



where we have neglected the fast oscillating terms that contribute to £ e g. For the values quoted in Figure l(a 
S eS = (15/100)7. 

The second method for coherent transfer applies the ideas of dark states in electromagnetically induced transparency 
(EIT) [TlJ to the collective states of two atoms. This method of transfer is remarkably robust, and allows for high 
fidelity rotations at large interatomic separations. If the incident fields are set to resonance (us — 0) and the fast 
oscillating terms in Eq. (fT3|) are ignored, then 

M = £„|c>-£» (16) 

is an (unnormalised) eigenstate of H e e that satisfies -ff ff = 0. The terms proportional to A can be ignored for 
approximately A > 2£ v , 2£^. Similar terms were ignored (with numerical justification) in Ref. [lfj]. For slowly- varying 
fields, Eq. ITS"]) can be written 

\tp(t)) = cos9(t) |c) - sin0(i) |b) , (17) 



where tan 6(t) = £ fJ- (t)/£ lJ (t) and slowly-varying means 9 -C <J £^(t) + £„(t). State \ift(t)) allows for pulsed population 
transfer between |c) and |b). We solve the full Hamiltonian numerically, and find that \tp(t)} is a good approximation 
to an eigenstate of_ff c ff. The population of levels |b) and |c) for the canonical counter- intuitive pulse sequence is 



shown in Figure 1(b) As noted in the caption, the transfer was possible at separations ~ (1/5)A . This regime is 
accessible to present experiment [2l|. Note the time taken for rotation is large, because when £ = 1.25, A ~ 1.2-y. 
In order to coherently transfer population between the Bell-states, the Rabi frequencies of the incident fields are 



modelled as sums of exponentials [Figure 2(a)| . The population of all four levels is shown in Figure [2(b) | which shows 



that high-fidelity transfer is possible using time-varying incident fields. The time taken for rotation is comparable 
with the emission time from a 2LA in free space. 



B. Emission properties 



We are interested in the effect of incident fields on the collective spontaneous-emission time and pattern from two 
atoms. We focus on the emission pattern from the single-excitation subspace, whose state is given by 



m)) = ^ iHdtt i-0(o)> , 



(18) 
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where H c g = Hs + Hi , \iIj(0)) is a linear combination of |b) and |c), and is unnormalised. The angular distribution 
of the emitted photon is |22| 

POO 

$(6>) sin 6d9 = / &t{{[){t)\&(0)S{e) |-0(*)> . (19) 



We denote the angular distribution affected by an incident field with <j>^(#). Without an incident field, the angular 
distribution of emission from states |b) and |c) can be written [22] 

$ b (6>) = L>(6»)sin(7rscos6») 2 , (20) 
$ c (6») = D{6) cos(tts cos Of. (21) 

Figure [3] shows the angular distribution of emission from |b) and |c) for an interatomic separation £ = 1/5. We have 
normalised the distributions so that $b(#) + 3>c(#) = D(8). 

In light of the manipulation methods proposed in Sec. IIII A1 we propose the ansatz 

Hi = £1 |c) (b| + H.c, (22) 

where fl is the Rabi frequency. Figure [3] compares the angular distribution when using the ansatz to the exact 
numerical solution. For the parameters quoted in Figure [3j the ansatz is a good approximation. The Rabi frequency 



Q = tanhf in order to account for the dark state preparation time shown explicitly for small t in Figure 2(b) 



The angular distribution with an incident field is not the dipolc emission pattern. The field causes the emission 
patterns from |b) and |c) to mix. Using the ansatz and 

\ip(t)} = cosOi|c) -isinfK|b), (23) 

the emission pattern for \ip(0)) = |c) is 

$o (6) sin 9d9 = J™ dt m ^+^f 7t [Sin cos | cosh - (T cos C - + sin |) sinh ^f] 2 , (24) 

for = Vr 2 — and ( = 2irscos9, and where we have normalised so that when O — s- 0, Eq. (f2Tj) is recovered. 
The matrix element 

(ip(t)\S\9)S(d)\^(t)) = 2~f D (9) cos(nt + its cos d) 2 , (25) 

for \ip(0)} = |c), explicitly shows the dependence of emission direction on il. The average of Eq. ([23)) for all times 
gives the (unnormalised) angular distribution of emission. Note that the emission is not simply $b|cb(0| 2 + < 5 ) c|c c (i)| 2 
for Mf)>=cb(t)|b} + Cc(t)|c}. 

The emission time is also altered by the presence of an incident field. Explicitly including decay, for Q 3> T, 

|V>(*)} =e~* 7/2 (c c cosm- ic b sinf2t) |c) + e~* 7/2 (c b cos fit - ic c sin Qt) |b) , (26) 

where IV'(O)) = Cb |b) +c c |c). This shows that for large O, the decay rate of two atoms becomes equal to that expected 
from a single atom. 

We have shown that for two atoms the emission direction and time is altered by an incident field that rotates within 
the one-excitation subspacc. Due to the speed and fidelity of the rotation b-c using pulscd-population-transfcr, and 
the large separations at which it is possible, we hope these effects are observable with present experiments. 



IV. THREE ATOMS 



The spontaneous emission from two atoms in the presence of coherent rotations allows some insight into the effect 
of incident fields on emission characteristics. However, two atoms is too few for both observation of directional 
superradiance and for exploiting collective decay properties for quantum information processing. So, we extend the 
analysis to three atoms. 

The three atoms are arranged in a linear configuration according to t*i = 0, = r and r$ = 2r. The energy 
level scheme is shown in Figure |4] Note that the definitions of |a) , |b) , |c), and |d) have altered from Sec. IIII1 The 
symmetric states are labelled |d) and |g). For quantum information processing, a qubit is most sensibly encoded using 
{|b) , |c)} as logical states. These are the slowest decaying excited states in the eight-level system, and in the limit 
£ — > have infinite decay times. They form the lower states of a decoherence-free subsystem Q. 
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FIG. 4: Energy level scheme of three atoms, with A± = (l/2)(<5 ± 3Ai 3 ), for 5 = V 8A ?2 + Af 3 . In the Dicke limit [H], 
|b) = (l/v / 6)(-2|001) + |010) + 1 100)) and |c) = (1/V2)(|010) - |100>). 



A. Manipulation 

Recently, there has been a proposal for the preparation, rotation, and readout of a qubit encoded across {|b) , |c)} 0|. 
So, wc assume that |b) can be prepared to high fidelity, with the effective coupling in the interaction picture with 
respect to the Hcrmitian part of H§ given by 



=e -iV^/i_^i|i(3, 



for £ M = d-E^, 5 = ^/8A^ 2 + A^ 3 , and K = 2A^ 2 +Ai 3 (Ai 3 — S). The laser frequency is w p = ui^ = (l/2)(Ai3— S)+cjq. 
Population transfer can be performed by using resonant fields and counter-intuitive pulse sequences. For this method, 
the incident helds resonantly couple b-g and c-g. So, by choosing the laser frequencies io v = ujt, g = 5 + coq and 
= ^c g = (l/2)(3Ai3 + S) + oj n in the interaction picture with respect to the Hermitian part of Hs the coupling c-g 



£ = Ai_ 2 (e? + 3A 13 )V<5 + A 13 r- i{ik . r+3Al3t)/2 , 2zk . r _ l){e m,2 g e 3iAl3i / 2 £ M )l , (28) 
cs 2V27[2A2 2 + A 13 (,5 + A 13 )] L V A ^' K ' 

for k v = = k, and the coupling b-g is 

£ * - «^iTA5i H 2A 'i2 + A? 3 )> - ^ + .cos, • r], (29) 

for K = 2A^ 2 + Ai 3 (Ai 3 — 8), and ? = 2Ai2Ai 3 (6A^ 2 — <5 2 + 3Aj 3 ). These couplings are such that, with judicious 
choice of Rabi frequency, they permit coherent Rabi oscillations between b-g and c-g. Setting all other couplings to 
zero, the dark state is written 

= £ cg \b) - S bg \c) . (30) 
As for the two-atom case, this is not an eigenstate of the total held Hamiltonian. Figure [5 (b)| shows the population of 



all eight levels when the pulse sequence shown in Figure 5(a) is applied to the collective system. Note that for three 



atoms, the laser helds are co-propagating so r |j k. Surprisingly, even in the presence of all the other couplings, there 
is a regime that admits coherent population transfer using counter-intuitive bichromatic fields. 

Similar to Sec. IIII A1 the manipulation can also be performed using collective two-photon Raman transitions 0]. 
Both methods of rotation rely on the splitting due to the dipolc-dipole interaction being large enough to enable 
incident fields to isolate the collective transitions. So, for experimental implementations, closely-spaced 2LAs are 
required. 
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FIG. 5: (a) Rabi frequencies of £ M = E n =i, 2.5, 3.25,4.8,. 207exp[-(9/100)(t - nt„) 2 /T 2 ] and £„ = 

E„=i, 1. 55,2.65... 207exp[-(9/100)(£ - nU) 2 /T 2 ] for T = 1.157" 1 , i M = 10T, U = i P + 3.75T, and laser frequencies of 
= tJbg = S + loq and lj m = uj cg = (l/2)(3Ai3 + S) + luq, and (b) populations of levels a, . . . , h when pulses of type (a) are 
applied to the collective eight-level system that is initially in |b) for a = ir/2 and the separation between neighbouring atoms 
is I = 1/5. 



B. Emission properties 

Using the ansatz 

Hi = n |c) (b| + H.c, (31) 

where il is the Rabi frequency, the decay of the levels |b) and |c) for 7b, 7c is cxp[— (l/2)(7b + Jc)t], where 7b 
and 7 C are the decay widths of |b) and |c) respectively. 

Comparing S(9) \ip(t)}, where \ip(t)} = cxp[-iH cS t}\ip(0)) for H cB = Rj from Eq. (31]) and |V>(0)) = |b), with 
S(9) |b) allows us to quantify the effect of the rotation on the emission direction. So, 

S(9) |b) =e^V^h^_^3 (K _ eKv) |a) (32) 




05 1 ~1.5 ~2 2l 3 

e 



FIG. 6: Average angular distribution of emitted photon with and without an incident held for a = 7r/2, £ = 1/5, and 
< t < 70y -1 , for emissions every ^St = 0.01, and for | -0(0)) = |b), normalised to the emission rate. 
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where r\ = (A i2 + A 13 )(<5 - 2A i2 - Ai 3 ), and 



m W)) =e 




(k - e K rj) cos nt + iV2(e K - l)«;sinfii ) |a) 



(33) 



The emission pattern has been altered by the field Hamiltonian, and includes a contribution from 



S(B)\c) 



(34) 



dependent on the time of emission and the size of f2. Figure [6] shows the numerical solution for the emission direction 
for the incident fields shown in Figure [5] The directional emission is proportional to the populations in |b) and |c). 
So, the directional emission is affected by the incident fields similarly to two 2LAs, and thus serves as a measure of 
the efficacy of population transfer. 



We have applied the ideas of STIRAP to multi-atom systems, and found that counter-intuitive pulse sequences can 
be used to transfer population, even if there is no null eigenstate of the field Hamiltonian. For two 2LAs, this method 
of transfer is robust at large separations, and so particularly amenable to experiment. 

For two and three 2LAs, we showed how the direction and time of the spontaneous emission is altered in the presence 
of an incident field. We gave analytical expressions for the directional matrix elements that depend on the strength of 
the incident fields. For three 2LAs we showed how the directional emission can give an indication of the efficacy of a 
rotation within a decoherence-free subsystem. We hope the results presented here will be useful for proof-of-principle 
experiments involving small numbers of dipole-coupled 2LAs. 
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by Macquarie University. 



[1] M. A. Rowe, D. Kielpinski, V. Meyer, C. A. Sackett, W. M. Itano, C. Munroe, and D. J. Wineland, Nature 409, 791 
(2001). 

[2] L. F. Wei, S. Y. Liu, and X. L. Lei, Phys. Rev. A 65, 062316 (2002). 

[3] C. Monroe, D. M. Meekhof, B. E. King, W. M. Itano, and D. J. Wineland, Phys. Rev. Lett. 75, 4714 (1995). 
[4] A. Steane, Appl. Phys. B 64, 623 (1997). 

[5] R. G. DeVoe and R. G. Brewer, Phys. Rev. Lett. 76, 2049 (1996). 
[6] V. Gorelik, K. Maeda, and H. Murai, RIKEN Review 44, 97 (2002). 
[7] Z. Ficek and R. Tanas, Phys. Rep. 372, 369 (2002). 

[8] J. Kempe, D. Bacon, D. A. Lidar, and K. B. Whaley, Phys. Rev. A 63, 042307 (2001). 

[9] P. G. Brooke, Phys. Rev. A 75, 022320 (2007). 
[10] A. Beige, S. F. Huelga, P. L. Knight, M. B. Plenio, and R. C. Thompson, J. Mod. Opt. 47, 401 (2000). 
[11] K. Bergmann, H. Theuer, and B. W. Shore, Rev. Mod. Phys. 70, 1003 (1998). 

[12] H. J. Carmichael, An Open Systems Approach to Quantum Optics ( Springer- Verlag, Berlin, Germany, 1993). 
[13] H. J. Carmichael and K. Kim, Opt. Commun. 179, 417 (2000). 

[14] A. A. Belavkin, B. Y. Zeldovich, A. M. Perelomov, and V. S. Popov, Sov. Phys. JETP 56, 264 (1969). 
[15] R. H. Lehmberg, Phys. Rev. A 2, 883 (1970). 
[16] G. S. Agarwal, Phys. Rev. A 2, 2038 (1970). 

[17] Q. A. Turchette, C. S. Wood, B. E. King, C. J. Myatt, D. Leibfried, W. M. Rano, C. Monroe, and D. J. Wineland, Phys. 
Rev. Lett 81, 3631 (1998). 

[18] C. F. Roos, G. P. T. Lancaster, M. Riebe, H. Haffner, W. Hansel, S. Guide, C. Becher, J. Eschner, F. Schmidt-Kaler, and 

R. Blatt, Phys. Rev. Lett 92, 220402 (2004). 
[19] Z. Ficek and R. Tanas, J. Mod. Opt. 50, 2765 (2003). 

[20] R. Tanas and Z. Ficek, J. Opt. B: Quantum Semiclass. Opt. 6, S90 (2004). 



V. 



CONCLUSION 



Acknowledgements 



10 



[21] J. Meijer, T. Vogel, B. Burchard, I. Rangelow, L. Bischoff, J. Wrachtrup, M. Domhan, F. Jelezko, W. Schnitzler, S. A. 

Schulz, et al., Appl. Phys. Lett. 83, 321 (2006). 
[22] J. P. Clemens, L. Horvath, B. C. Sanders, and H. J. Carmichael, Phys. Rev. A 68, 023809 (2003). 
[23] R. H. Dicke, Phys. Rev. 93, 99 (1954). 



